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Abstract: At energies much less than the electron mass m the effects of quantum fluctuations 
in the vacuum due to virtual electron loops can be included by extending the Maxwell Lagrangian 
by additional non-renormalizable terms corresponding to the Uehling and Euler-Heisenberg in- 
teractions. This effective field theory is used to calculate the properties of the QED vacuum at 
temperatures T <C m. By a redefinition of the electromagnetic field, the Uehling term is shown 
not to contribute. The Stefan-Boltzmann energy density is thus found to be modified by a term 
proportional with T 8 /m in agreement with the semi-classical result of Barton. The speed of light 
in blackbody radiation is smaller than one. Similarly, the correction to the energy density of the 
vacuum between two metallic parallel plates diverges like l/m 4 z 8 at a distance from one of the 
plates z — > 0. While the integral of the regularized energy density is thus divergent, the regularized 
integral is finite and corresponds to a correction to the Casimir force which varies with the sepa- 
ration L between the plates as 1/m L . This result is in seemingly disagreement with a previous 
result for the radiative correction to the Casimir force which gives a correction varying like 1 /mL 5 
f~| ■ in a calculation using full QED. 
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The use of effective field theories is steadily increasing in the investigation of quantum 
phenomena at high energies [Q] ||] ||] Q . Since these theories are generally non-renormalizable, 
^ it has taken some time before one has fully understood how they are supposed to be used 

in higher orders of perturbation theory. Their applications are only valid below a given 
energy scale. Loop integrations are usually divergent and bring in high-energy degrees of 
freedom which are not described by the effective theory. This problem is cured by regu- 
larization. It is most directly done by introducing an explicit momentum cutoff. Its effect 
can later be removed by the introduction of higher dimensional counterterms to cancel the 
divergences and yield finite results. An alternative and much more compact method is to 
use dimensional regularization. It seems to be ideally suited for effective field theories. 

At energies of the order of the electron mass m, the most important interactions are 
described by QED. In principle, one can also use this theory at much lower energies to 
explain phenomena in atomic physics or condensed matter physics. But that would be 
technically very difficult and one relies instead on other, effective theories involving more 
appropriate degrees of freedom. The effects of high energy physics are then coded into 
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the values of the coupling constants of new interactions. A typical example is NRQED||] 
which has been developed during the last few years |]]. It makes it possible to calculate 
with the highest accuracy radiative and relativistic corrections to the energy levels of the 
simplest atoms such as positronium[|jJ and muonium[||. 

In the QED vacuum there are no physical electrons or photons. Since the photons are 
massless they will dominate all virtual processes for energies much less than the electron 
mass. These quantum processes can therefore be described by an effective theory for 
the electromagnetic field where the effects of virtual electron loops appear as new, non- 
renormalizable interactions. The dominant terms are the Uehling corrections to the photon 



propagator || and the Euler-Heisenberg four-photon interaction [10]. This effective theory 
is discussed in the next section. We will show that the Uehling interactions have no effects 
in the QED vacuum and can be transformed away by a redefinition of the photon field. 
This simplification is further discussed in Appendix A. 

When the temperature is non-zero and below the electron mass m, we can use the 
effective theory to calculate the free energy of the photon field. This is done in Section 2 
where first the radiative correction to the free energy of interacting photons is derived. 
Using thermodynamics we then also have the lowest correction to the Stefan-Boltzmann 
law for black body radiation. It results from the Euler-Heisenberg interaction and is found 



to be in agreement with a previous, semi-classical calculation by Barton] 11]. In full QED 
this correction is due to very complicated 3-loop diagrams. On the other hand, in the 
effective theory it is directly given by a product of two 1-loop diagrams which are straight- 
forward to evaluate. This clearly demonstrates how powerful effective field theories can 
be for calculational purposes. The energy-momentum tensor for the interacting system 
is constructed and its expectation value is found to be in agreement with the thermody- 
namic derivation. A few technical aspects of this calculation can be found in Appendix 
B. Also the temperature-dependent corrections to the electric permittivity and magnetic 
susceptibility of the vacuum are calculated from the electromagnetic polarization tensor. 
It gives a reduction of the speed of light in blackbody radiation in agreement with a recent 



result by LaTorre, Pascual and Tarrach[12 



Here we consider only photons at low temperatures. In the opposite limit, i.e. T S> 
m, one can also use the methods of effective field theory to simplify the calculations of 



thermodynamic properties] 13 1. The characteristic scale is then the temperature and the 
electron mass can be set equal to zero. For length scales > 1/T it is then safe to neglect 
the modes with non-zero Matsubara index and the effective theory becomes 3-dimensional. 
The coefficient of the Stefan-Boltzmann energy gets then radiative corrections and there 
will be additional terms in the free energy which now varies logarithmically with the 
temperature p 



The Casimir effect for two parallel plates is considered in Section 3. After having con- 
structed the electromagnetic multipoles in the Coulomb gauge, the free field is quantized. 
Quantitative results for the field fluctuations are then given by the different correlators. 
In Appendix B they are also calculated from the propagator of the photon field confined 
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between two parallel plates as constructed by Bordag, Robaschik and Wieczorekj^] in 
the Lorentz gauge. The correlators are in general found to be strongly dependent on 
the position between the plates. However, the vacuum energy density in the free theory 



is constant and gives the well-known Casimir force[17| [18| |19|| . It has recently |20|| been 
measured with much greater precision than previously [21]. However, when we include 
the Euler-Heisenberg interaction in the effective Lagrangian, we find that the regularized 
energy density diverges near the plates. It thus results in an infinite Casimir energy when 
integrated over the volume between the plates. This physical meaningless result can be 
avoided by first integrating and then regularizing the Casimir energy. We then obtain a 
finite result for the correction to the Casimir force which is found to vary with the dis- 
tance L between the plates as 1/m L 8 . This is not in agreement with a corresponding 



calculation using full QED by Bordag et al\15\ who obtained a result varying like 1/mL 5 . 



Later this was confirmed by Robaschik, Scharnhorst and Wieczorek[16] using a similar 
approach. If the Uehling term had not have been transformed away, it would on the other 
hand have given a correction going like l/m 2 L 6 for dimensional reasons. In general one 
should obtain the same results from both the full and the effective theory as long as they 
are properly matched. 

It was first noted by Deutsch and Candelasp2j| that the operations of integration and 
regularization do not in general commute in this type of problems. In order to analyze 
this conflict more easily, a simpler model of an interacting, massless scalar field in 1+1 
dimensions has been investigated |2^j . The model corresponds closely to the above de- 
scription of interacting photons between two plates but without the uninteresting degrees 
of freedom in the transverse directions. It is then easy to show that the lack of com- 
mutativity between integration and regularization is due to the breakdown of different 
regularization schemes when one gets sufficiently close to one of the walls. In this region 
unknown, cutoff-dependent effects show up reflecting the microscopic properties of the 
physical walls. 



1 Interacting photons at low energies 

QED describes the interactions between the electromagnetic field A^(x) and the electron 
field ip(x) by the standard Lagrangian 

C(A, VO = -±F% + V47 M (^ - eAJ - rr# (1) 

where F^ v = d^A v — d u A^ is the field strength. At energies much below the electron 
mass m and with no matter present, the first term which is the Maxwell Lagrangian, will 
dominate and describe free photons in the quantized theory. The effects of the electrons 
are only felt through the coupling of photons to the virtual electron loops in the vacuum. 
These interactions are suppressed by inverse powers of the electron mass compared with 
the Maxwell term and are thus very small. 
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One can include the effects of these virtual processes by extending the Maxwell La- 
grangian to include higher dimensional operators that describe these additional interac- 
tions. Lorentz invariance combined with parity and gauge invariance require these interac- 
tions to be constructed from the electromagnetic field tensor F^ u and derivatives thereof. 
The coefficients of such operators in this effective Lagrangian C e ff can then be obtained 
by matching results with those obtained from full QED||. In principle this will lead to 
the same result as directly integrating out the electron field in the vacuum-to-vacuum 
transition amplitude or partition function, 

e iJd 4 xC eff (A) = fx>^Vi)e^ dAxC ^ A) (2) 



Since the QED Lagrangian is quadratic in the electron field, the integration gives a func- 
tional determinant and the effective Maxwell Lagrangian can formally be written as 

C eff (A) = -If% - .Tr log [VCify - eA,) - m] (3) 

It can expanded in an infinite series of operators O n with increasing dimensions D n > 4 
preceeded by expansion coefficients proportional with m 4 ~ Dn . The lowest order contri- 
bution is due the to 1-loop correction of the photon propagator as shown in Fig.||. It 
modifies the Coulomb potential at short distances. This is the Uehling interaction^] and 
it is quadratic in the photon field. Expanding it to lowest order in the inverse electron 
mass it gives rise to the interaction p4| 

a 

607T771 2 



Lv = T^F^F^ (4) 



where a = e 2 /47r is the fine structure constant and □ = dnd^. Inserting extra photon 




□ D 2 



Figure 1: Uehling corrections to the photon propagator. 

lines into the electron loop gives 0(a) or higher corrections to the coefficients of these 
interactions which will be neglected in the following. Expanding to higher orders in 1/m, 
one finds in addition to the dimension D = 8 Uehling term, the Euler-Heisenberg effective 



interaction 1 1C] 



C-eh 



a 2 
90m 4 



(5) 



where F^ = 2 e ^p^F pcr is the dual field strength. It can be derived analytically from the 
functional determinant in (|3[) when the field is constant as shown in many textbooks p4|1. 
This effective interaction gives a point coupling between four photons. In full QED it 
arises from the coupling of the photons to a virtual electron loop as shown in Fig|| To 
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Figure 2: The Eulcr-Heisenberg four- photon vertex. 



this order in the inverse electron mass we thus have the effective Lagrangian 



f 



C eff = --F^ + C u + C EH (6) 

for the interacting Maxwell field. The equation of motion for the free field is = 0. 

It can now be used to simplify the interaction terms which in our case with no matter 
means that the Uehling term Cjj can be effectively set equal to zero. Since the equation 
of motion is only satisfied by on-shell photons, one may then question the validity of 
this simplification where the interactions are used to generate loop diagrams with virtual 
photons. A justification for this procedure is given in Appendix A where we for simplicity 
consider a similar scalar theory which also couples to matter. Here it can be said that 
the effective Lagrangian is supposed to be used in the functional integral for the partition 
function where one is generally allowed to shift integration variables by local redefinitions 
of the field. In fact, we see that by the transformation 

A u ^A u -\ ^-—UA U (7) 

M M 307rm 2 " V ; 

the Uehling interaction is eliminated. The field redefinition does not modify the four- 
photon (||) interaction to this order in the electron mass. 

The effective theory of photons at low energies based on the Euler-Heisenberg La- 
grangian is a non-renormalizable theory since the interaction has dimension D = 8. Until 



recently it has therefore been treated as a classical theory. It seems that Halter [25] was 
the first to consider it instead as a fullfledged quantum theory in which one consistently 
can calculate radiative corrections. He considered in particular higher order corrections to 
photon-photon scattering and calculated the contribution from one-loop diagrams where 
the Euler-Heisenberg vertex appeared twice, i.e. giving an amplitude going like uj 8 /m s 
where uj is the photon energy. This corresponds to the single action of a dimension D = 12 
vertex and represents a subleading contribution^]. More important are the 0(a) cor- 
rections to the tree-level result due to higher order corrections to the coefficient of the 
Euler-Heisenberg interaction (||). These where originally calculated by Ritus[27| and later 



confirmed by Reuter, Schmidt and Shubert[28|. Secondly, as pointed out recently by 
Dicus, Kao and Repko]^s|l, the dimension D = 10 interaction in the expansion of the ef- 
fective theory will become important at sufficiently high energies and thus also contribute 
at tree-level. This new interaction, which they have constructed, involves four photon 
field strengths with two derivatives. In the following we will ignore these higher order 
corrections. 
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2 QED at low temperatures 



The energy density of free photons at non-zero temperature T is given by the Stefan- 
Boltzmann law 

^2 

(8) 



7T" 



£ = — T 4 
15 



Due to interactions with virtual electron-positron pairs this fundamental formula will in 
general have quantum corrections. The lowest contribution results from QED diagram in 
Fig.^ which can be obtained from the photon self energy in Fig.[l] by connecting the two 
photon lines. It is therefore the Uehling contribution to the vacuum energy. However, 







Figure 3: Possible Uehling correction to the free energy to order 1/m 4 . 

when it is calculated for temperatures T <C m it is found to be exponentially suppressed 
like exp (— m/T) [pC|1 . This is consistent with our observation that the Uehling interactions 
can be transformed away in the low-energy effective theory. 

The first correction to the Stefan-Boltzmann law is therefore given by Euler-Heisenberg 
interaction and must for dimensional reasons vary with the temperature like T 8 /m 4 . We 
will calculate its absolute magnitude using the imaginary-time reformulation of the effec- 
tive theory. In addition to an over-all change of sign, the part involving the dual tensor 
in the effective Lagrangian will also change sign so that the we have the Euclidean theory 

„2 



C E 



a 



90m 4 



-(E 2 + B 2 ) 

2 V ' 



2a 2 
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45m 4 



E 2 + B 2 ) 2 -7(E-B) 2 



(9) 
(10) 



For the purpose of doing perturbation theory we need the correlator (F^ u (x)F a p(y)) of 
the free field. It follows from the free photon propagator with 4-momentum k^ = (k, k^) 
in Euclidean space, 



(A„(k)A v (-k)) = ± 



<W + (£ - 1) 



kuk v 



k 2 



where we can choose the gauge- fixing parameter £ = 1. Then we have 

i 

{A^(k)F a p(-k) ) = ~2 [k a 6fj,p - kp5^ a ] 



(11) 



(12) 



which will be needed when we calculate the polarization tensor. The gauge invariant 
correlator of the field tensor is therefore 



(F lxu (k)F a/3 (-k)) = 7£ [k^kpSv 



k^ k a 5 v p 



(13) 
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When the photon field is in thermal equilibrium at temperature T, the fourth component 
of the four-momentum vector fc„ is quantized, k± = u) n = lixTn where the Matsubara 
index n takes all positive and negative integer values for bosons. Loop integrations are 
then done by the sum-integral 

00 r H 3 k 

T Z W? (14> 

n=— 00 v / 

where the 3-dimensional momentum integration is dimensionally regularized and the Mat- 
subara summation is regularized using zeta-functions. Many such sum-integrals have been 
calculated by Arnold and Zhaij3l| and Braaten and Nieto|l3|]. Some of the needed integrals 
are given in Appendix B. We notice here that integrals on the form JJ. k 2n are non-zero 
only when n < 0. When n > they look like being highly divergent, but are actually zero 
when using dimensional regularization. 

As a check, we can now obtain the energy density of free photons which in the Eu- 
clidean formulation is given by 

£(x) = -i(E 2 (x)-B 2 (x)> (15) 

Using the basic integral 

§£^ = ^(V- 4 W*) (16) 

we obtain from (12) above 

4 

and 

2k 2 7T 2 7 ' <l 

v 4 

which reproduces the Stefan-Boltzmann energy (|8|). The corresponding free energy density 
is 



T=-^—T A (19) 
45 

and equals the negative of the pressure P. 
2.1 Corrections to the free energy 

The free energy of interacting photons follows directly from the partition function 

Z = e -P!**T = fvApe- ^ dix Ce (20) 
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where f^d^x = Jq dx\ J d?x. To lowest order in T/m the correction to the free energy 
density is thus simply given by the expectation value AT = (ACe) in the free theory 
where ACe is the Euclidean Euler-Heisenberg interaction (g). For the evaluation of the 
expectation value, it is convenient to rewrite it as 

a 2 5 
AT = QQ m 4 ( fFfwFvpFpaFam — -F^ u F u ^Fp a F a p ) (21) 

after having expressed the dual field tensor in terms of the usual one. We see that the 
correction to the free energy corresponds to the calculation of the 2-loop Feynman diagram 
in Fig.0. It is seen to be really the product of two 1-loop diagrams and therefore almost 




Figure 4: Euler-Heisenberg correction to the free energy. 

trivial to calculate compared with the corresponding 3-loop diagram in full QED. Using 
the correlator ( |l3| ) and doing all the momentum contractions, we find the result to be 
given by the double sum-integral 



AT 



22a 2 rr(p. q ) 



45m 4 T P Tq p*q 

which can be evaluated from the basic integral ( |i~6|) which gives 

p 2 q 2 675 

The free energy density of interacting photons to lowest order in perturbation theory is 
therefore 

^_ 7T 2 4 227T 4 a 2 r 8 

45 3 5 -5 3 m 4 ( ' 

which is also the negative pressure. 

This result has previously been derived by Barton [JO]] using a semi-classical method 
and treating the interacting photon g material medium. Since the entropy is given 

by the derivative of the free energy with respect to the temperature, the energy density 
follows from £ = (1 - Td/dT)T as 

vr 2 „ 7-22 ?r 4 a 2 T 8 
£ = — T 4 + - f " — A (25) 
15 3 5 • 5 3 m 4 y ' 

Obviously the corrections due to the Euler-Heisenberg interaction are totally negligible at 
ordinary temperatures. One can perhaps speculate that they may be of relevance under 
special astrophysical conditions. In this spirit Barton fi~l"[| has investigated the implications 
of these photon interactions for the Planck distribution of blackbody radiation. 
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2.2 The energy-momentum tensor 



The previous calculation was essentially based upon the thermodynamics of the interacting 
photon gas. In particular we see that the factor seven which relates the correction to the 
energy to the free energy comes from taking 1 — Td/dT acting on T 8 which seems to be 
of purely thermodynamic origin. However, we also know that the energy density £ and 
pressure P are both directly given by the energy-momentum tensor for the interacting 
system. It is of some interest to see how this field-theoretic approach can reproduce the 
above thermodynamic results. 

Placing the interacting theory (||) in a curved spacetime with metric g^ u , the energy- 
momentum tensor is in general given by 

nn , s 2 8S[A] , , 

T u Jx) = -z — (26) 

where S[A] is the corresponding action functional for the system and g = det(g fiu ). The 
result in flat spacetime is given by the standard tensor 

Tj£ = F fl \F u \ — -S^F aj3 F a j3 (27) 
in free Maxwell theory plus a correction 

„2 



1 . 
+ 2 ^ 



(F a pF a(3 ) 2 - 7 -{F aP F a(] ) 2 



(28) 



due to the Euler-Heisenberg interaction. 

In the following it will be convenient to express the energy-momentum tensor T^ u = 
T*5 + T^av i n terms of the components of the electric and magnetic fields. In particular, 
for the calculation of the energy density £ = ( T44 }j„t where the expectation value is taken 
in the interacting theory, we need the Maxwell result 



\ (E 2 - B 2 ) (29) 

and the Euler-Heisenberg contribution 



rpM 

U4 2 



T; 



EH 
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2a 2 
45m 4 



(B 4 - 3E 4 + 7(E • B) 2 - 2E 2 B 2 ) (30) 



Using the previous free equal-time correlators (EiEj) = —(BiBj) = — <5jj7r 2 T 4 /45 and 
( EiBj ) = combined with Wick's theorem, we obtain the correction 

,^ PH , 22vr 4 a 2 r 8 . , 
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Figure 5: Feynman diagram giving the expectation value in the free theory of the interacting 
energy-momentum tensor 



due to the photon interactions. The result corresponds to the evaluation of the Feynman 
diagram in Figj^ which in the previous section appeared in the calculation of the correc- 
tion to the free energy. But there is an additional contribution coming from the photon 
interactions as is seen from the general expression for the expectation value 



iiit 



| J VA^ (T44 + Tg H ) e~ ^ d " xCl 
Jo 



(Tj£) + (T£») + A(T i 



EH < 



iM • 
44 . 



(32) 

(33) 
(34) 



when expanded to lowest order in perturbation theory. Writing out the last term, it is 
seen to be given by 



a(th: 



2a 2 



f[E 2 (y) + B 2 ( 2/ )] 2 -7(E(y).B(y)) 2 



(35) 



It is a non-local contribution receiving contributions from all over spacetime and corre- 
sponds to the Feynman diagram in Fig.|6] with three propagators. This makes now the 
calculation more difficult than in the previous case of the 2-loop diagram in Fig.||. With 




Figure 6: Feynman diagram for the non-local contribution to the expectation value of the energy- 
momentum tensor 

most of the needed details of the calculation in Appendix B, we find the result 



A(T 44 ) 



n A a 2 T 8 
3 4 • 5 3 m 4 



10 



10 

y 



50 „ 

Y + 7 



+ 2 -4 + (-10) -7 - U^j 

6-22vr 4 a 2 T 8 



3 5 -5 3 



(36) 



Adding this contribution to (|3l| ) we obtain the previous thermodynamic result fl25J) for 
the correction A£ to the Stefan-Boltzmann result. 



While the energy-momentum tensor (|27| ) for the free Maxwell field is traceless, im- 
plying the relation £ = 3P between energy density and pressure, the Euler-Heisenberg 
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interaction ACe gives a contribution (|2g) which has a non-zero trace. Actually, we now 
find 

= AAC E (37) 

Since the correction to the pressure due to the interaction is AP = —( ACe )> we then have 
the relation £ — 3P = 4AP which is just the previous thermodynamic result AS = TAP. 



2.3 Response functions 



The polarization of the vacuum is described by the photon selfenergy IL^A;). In lowest 
order perturbation theory it is given by the correlator ( A^{k)ACE{^)A u {— k) ) correspond- 
ing to the 1-loop Feynman diagram in Fig.^. At zero temperature it has been calculated 




Figure 7: Euler-Heisenberg correction to the photon propagator. 



by Halter pq| who showed that the photon remains massless as expected. Here it is corre- 
spondingly given by the sum-integral 



tea 1 



1 r 



(A; • q)\ v + k 2 (k ■ q)(k fM q u + q^ky) 



45m 4 J q q 2 

All the terms can be evaluated from the basic sum-integral (|l6| ) which gives 

44vr 2 a 2 T 4 r 2 

= 2025 m 4 r ^ ~ 26 ^ 6u ^ 

— k^k v — 2k 2 6 hi/ + 2k4{k ll 5 U 4 + fe I/ <5^4 y 



(38) 



(39) 



Introducing the transverse projection tensors P^v{k) and Q^ u (k) as defined by Weldon 
||32j| we can write the above result as 



U fU/ (k)=U T (k)P /a , + U L (k)Q 



(40) 



in terms of the scalar response functions nj'(fc) and Yli,{k). From the explicit form of the 
projection tensors, we then have in general 



n L (A0 



k 2 



n 44 , u T {k) = -(u^ + u L ) 



From our result (39) we then get for these functions 

44vr 2 a 2 T 4 



n L (fc) 



2025 m 4 



(41) 



(42) 
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and 



The transition back to Minkowski space is now achieved by the standard analytical con- 
tinuation z/c4 — * uj for the time component of the momentum vector while the space 
components remain unchanged. In this way we obtain the physical response functions 
7Tt,l(cJ) k) from the above results. We then have the photon propagator in Minkowski 
space from which we can read off the frequency dependent permittivity e(oj, k) and sus- 
ceptibility k). One then finds p2[ 

e(u,k) = l-^ (44) 

where now k^ = (lo, k) with k 2 = u> 2 — k 2 is the physical 4-momentum in Minkowski space. 
Similarly, 

1 k 2 TT T -UJ 2 7r L 

1 + ( 45 ) 



fi(w,k) k 2 k 2 
In the static limit uj — > we thus obtain for the interacting photon gas 



44vr 2 a 2 T 4 , , 

a=l + j 46) 

P 2025 m 4 v ; 



so that the speed of light becomes 

_ 447r 2 a 2 T 4 




(47) 



2025 m 4 

It is reduced below its ordinary vacuum value due to the interaction of the photons with 



virtual electron-positron pairs as first calculated by Tarrach|33||. A corrected result has 



since then been derived in full QED from 2-loop diagrams[12] and agrees with the above 
result obtained more directly in the effective theory. 



3 The Casimir Effect 



In the previous section where the photons were at equilibrium at finite temperature, the 
corresponding electromagnetic field was periodic in the imaginary time dimension. It gives 
rise to a finite, temperature-dependent vacuum energy. This situation is mathematically 
identical to compactifying one of the spatial dimensions instead. For instance, one can 
imagine the field at zero temperature confined by boundary conditions in such a way that 
it is periodic in the z-direction with period L. The result for the free energy calculated 
in the previous section can then be used to obtain the vacuum energy for this confined 
field by the substitution T — * 1/L. In the case of free photons, the Stefan-Boltzmann 
free energy density (U) gives the vacuum energy density £ = — 7r 2 /45L 4 . This is now a 
quantum effect due to the vacuum fluctuations of the field. 
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3.1 Electromagnetic multipoles and quantization 



Vacuum energies induced by non-trivial boundary conditions in space imposed on a quan- 



tum field is now generically called Casimir energies and have been much studied [ 18 |[|T 



We will here consider the original Casimir effect[17] due to fluctuations of the electromag- 
netic field A^(x) = (Ao(x, t), A(x, t)) confined between two parallel plates with separation 
L and with normal n along the z-axis. We take the plates to be perfect conductors so that 
the fields satisfy the metallic boundary conditions 

nxE = n- B = (48) 

at all frequencies. For this problem we find it most convenient to work in the Coulomb 
gauge V • A = so that the Maxwell equation d^F^ = separates into the wave equation 

(d 2 -V 2 )A(x,t) = (49) 

for the magnetic potential and the Poisson equation V 2 ^4o( x , t) = for the electric po- 
tential. Since the latter is zero on both plates, we can then simply take Aq = 0. Thus we 
have the field strengths E = — A and B = V x A. 

The magnetic potential A(x, t) is most easily quantized by expanding it into electro- 
magnetic multipoles |3^| . From the symmetry of the problem we know that the solutions 
of the wave equation ( f49| ) with a definite frequency u> will be labeled by wave vectors 
k = (k.T,k) with lo = |k|. Here = (k x ,k y ) represents the transverse components and 
k = k z is the longitudinal component. Transverse electric (TE) modes will have E z = 
and can thus be written on the form 

A*(x) = V x u£(x) (50) 

while the transverse magnetic (TM) modes have B z = so that 

Af (x) = V x [V x Uf (x)] (51) 

Using the wave equation ( |4^ ) we have for the corresponding potentials 

U^(x) = Nnsm(kz)e ikT ^ T (52) 

Uf(x) = — ncos{kz)e ikT -* T (53) 
no 



where xy = (x,y). The boundary conditions ([48|) are now satisfied when taking 

k=yn, n = 0,1,2,... (54) 

The corresponding frequencies are therefore also quantized, 

u = (]4 + (nvr/L) 2 ) 1 / 2 (55) 
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We see that the special mode with n = has B z = and is thus a TM mode. In this way 
we have found all the multipoles except for their absolute size given by the constant N. 
It can be determined from the normalization condition which we take to be 



J d 3 x A^ kT (x) • A^(x) = 5 XX , <W (27r) 2 6(k T - k' T ) 
where A = (E,M). After some integrations, one obtains 

n = 



(56) 



N 



j_ l 

k T V L ' 
1 l~2 



L - 



n 



1.2, 



(57) 



It is the n = TM mode which has the exceptional normalization. The full vector potential 
can now be expanded in terms of these transverse modes as 



A(x,t) 



E E 

X=E,M ra=0 



d 2 k 7 
(2vr) 2 V 2w 



n A A A 



(x)e 



-itat 



+ C.C. 



(58) 



In the quantum theory the expansion coefficients a„ k are annihilation operators satisfying 
the canonical commutator 



A T ,4l T ] = S xy <W (2vr) 2 5(k T - k^) 



(59) 



Since the vacuum state by definition is annihilated by these operators, the expectation 
values of the electric and magnetic fields in this state are zero. But the fields have non-zero 
fluctuations which will be calculated. 



3.2 Casimir energy 



The most important physical effect of the vacuum fluctuations will be the attraction 

It is 



between the plates discovered by Casimir [17] and now experimentally verified |2C 
due to the zero-point energy of the vacuum state and is given by the divergent expression 



1 ffkr_ 

2 J (2vr)2 



k T + 2 



oo 

E 

71=1 



kf, 



(60) 



The first term is due to the TM n = mode while the sum includes equal contributions 
from the TE and the TM modes. Obviously, both the sum and the integrals are strongly di- 
vergent. Usually they are reg ulated by a UV cutoff pi], but we will find it more convenient 
to use dimensional regularization of the transverse momentum integrations combined with 
zeta-function regularization of the longitudinal momentum summations as in the previous 
section. We will therefore again use the standard integral 



I{d,p) 



d d k T 
(2tt) c 



k d ~p r(gr) 
(4tt)<*/ 2 r(f) 



(61) 
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When k = the result is zero so that the first term in ( p0[) gives zero. In the second term 
we need 1(2, —1) = —k 3 /6n so that the Casimir energy becomes 

1 ^ (-K \ 3 vr 2 



n=l 



B = -^2, I") =-720i? (62> 



when we use the value £(— 3) = 1/120 for the Riemann zeta- function. This is the energy 
per unit plate area so that the energy density is £ = E/L = — 7r 2 /720L 4 . The attractive 
Casimir force between the plates, is given by the derivative of the energy (|62|) and thus 
becomes F = -vr 2 /240L 4 . We notice that the above energy density follows from the 



Stefan-Boltzmann free energy density (19) by the substitution T — ► 1/2L. Forcing the 
field to satisfy metallic boundary conditions therefore corresponds to taking L to be half 
of the period 1/T at finite temperature. 



3.3 Field fluctuations 

In the above we have implicitly assumed that the energy density is constant between the 
two plates. But there is no deep reason for that to be the case. For instance, the Casimir 
energy density inside a spherical, metallic shell increases outward from the center and 



diverges when one approaches the shell [35]. In the present case with parallel plates the 
energy density is in fact constant. However, as soon as we include field interactions we 
will find that the density varies with the position between plates. 

The electric field E = (Ey, E z ) is given by the time derivative of the magnetic potential 
(|5§|). For the fluctuations in the transverse components we then find 

while the longitudinal fluctuations are given by 

1 ~ fd 2 k T ( k 2 \ 2/I N 

^^)) = iE/(2#(--^j cos2 ^) ( 64 ) 

Corresponding expressions for the fluctuations in the magnetic field are also easy to derive 
and can in fact be directly obtained from the relations ( B 2 ) = (B 2 ) = —(El) and 
( B 2 ) = —(E 2 ). Using the regularized integral (|6l|) which gives 1(2, 1) = —k/2ir, we find 
that the transverse fluctuations become 

7T 2 °° 

( E|(x) ) = £ n 3 (l - cos 20n) (65) 

n=l 

where 6 = ttz/L is the dimensionless distance between the plates. Again with zeta- function 
regularization the first sum is £(— 3) = 1/120 while the last sum becomes 

00 /l d \ 3 00 1 

^n 3 cos2#n = - -— ^ sin 20n = -F(0) (66) 
n=l \2d0J n=1 8 



15 



where the function 

Id 3 3 2 

2 dO 6 sm 4 sur 6* 

gives the position-dependence of the field fluctuations, 

Similarly, we obtain for the longitudinal fluctuation (64) 

( £ 2 (x) ) = _ ^ n 3 (l + cos 29n) (69) 



which then becomees 



n=l 



(^W»'*' = iSi(iS + F «) <70> 

after regularization. These results now also describe the fluctuations of the magnetic 
field using the above relations to the corresponding electric fluctuations. Combining the 
transverse and longitudinal components we thus have for the fluctuations of the full field 
strengths 

<*"«>"•— ^(s-w). <B , wr , ~w(s+w) (71) 

In addition we have ( EiBj ) = which also follows directly from time-reversal invariance. 



These results have previously been obtained by Liitken and Ravndal f34j| using a dif- 
ferent regularization scheme formed by letting the difference between the time coordinates 
of the two fields in the correlators be imaginary and slightly non-zero in magnitude. More 



recently, Bordag et al [15| have considered the same problem in the Lorentz gauge and 
obtained a closed expression for the photon propagator in the presence of the two plates. 
Letting the two coordinates in the propagator approach each other, one obtains expres- 
sions for the local fluctuations of the quantum fields. This is done in Appendix C. After 
regularization, now based upon a point-split in the spatial direction, we then again recover 
the above results. 



Because of the boundary condition fl48|) the transverse electric field is constrained to 
be zero on the walls. The fluctuations in the transverse components must therefore be 
zero for z = and z = L which we see is consistent with the general expression (|63|). 
However, the regularized fluctuations (|7l|) are seen to diverge when one approaches one 
of the plates. When z — > the leading terms are 

< E2 < X »=° = I^' (B 2 (x)a = -^ (72) 



z' 
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These expressions will represent the field close to any metallic surface, plane or curved, 
since when one gets close enough, it will always appear as plane. Since the transverse 
fluctuations are exactly zero on the plate and highly divergent just outside, it is clear that 
the physics very close to the plate is complicated. In fact, it will depend on the microscopic 
properties of the wall reflected in the existence of a non-zero cut-off in this region. In other 
words, when one gets very close to a wall, the detailed properties of the fluctuations will 
depend on this cut-off and the regularization becomes dependent on which scheme is being 
used. Using for example a non-zero point-split e in the imaginary time dimension, it is 
easy to show|2^] that the limit e — > cannot be taken when the distance to the wall is less 
than e. And the magnitude of e in our effective field theory approach is set by the inverse 
electron mass. 

The position-dependence of the correlators ( E 2 ) and ( B 2 ) can in principle be mea- 
sured. An atom outside a single plate will be attracted to the plate by the Casimir-Polder 



force [36] which has a position dependence directly given by the one in the correlators] 34]. 
This force, or its equivalent van der Waals force at shorter distances, has now been ex- 
perimentally verified [37]. The electric fluctuations (E 2 ) will affect the energy levels of an 



atom near a wall [ 34 1. By moving the atom around one can then map the spatial variation 
of the vacuum fluctuations by spectroscopic methods ]38[]. One will then observe that these 
effects become much stronger as one gets very near the plate or wall. But the obtained 
results for the local fluctuations are not valid when one gets too close to a wall. Then new 
physics determined by microscopic properties of the wall will start to become relevant and 
are outside the present description. 



3.4 Euler-Heisenberg interactions 



It is of interest to calculate the radiative corrections to the one- loop Casimir energy (62) 
due to the interactions of the fluctuating electromagnetic field with the virtual electrons 
in the vacuum. The first correction in QED comes from the Uehling interaction in Fig. 



which has been evaluated by Xue[39|. Instead of the physical boundary conditions (|4| 
he constrained the photon field to be periodic in the z-direction with period 2L. As 
previously noted, this gives the correct result for the 1-loop term. He finds that the 2- 
loop contribution is exponentially suppressed like exp (— 2mL) and therefore completely 
negligible in the physical limit L> \jm. Since the system with these boundary conditions 
is now mathematically equivalent to being at finite temperature T = 1/2L, this outcome 
is consistent with result of WoloshynpO| discussed in the first section. 



A much different calculation by Bordag et al and later by Robaschik et al [ |ilj| ] 
obtains a radiative correction to the Casimir force proportional with a/mL 5 using the 
physical boundary conditions (|4~8|) for the photon field. The electron field is assumed to 
be unaffected by the metallic plates. The correction is small and decreases faster with 
the plate separation than the leading term. But it is somewhat surprising to see that it 
depends on an odd power of the electron mass. Since this calculation is done within full 
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QED, it should be valid down to distances of the order of the inverse electron mass. So in 
this case one is no longer discussing a physical plate made up of atoms, but a more ideal, 
mathematical situation. It is thus not obvious that this result should agree with what 
is obtained using low-energy, effective field theory. Then the first radiative correction to 
the 1-loop result should result from the Uehling term in (Q). Since it is proportional to 
a/m 2 , it could give a correction to the Casimir effect varying like a/m?L e for dimensional 
reasons. However, since the Uehling interaction can be transformed away, the first non- 
zero contribution will result from the Euler-Heisenberg interaction (||). In Minkowski space 
it is given by the perturbation 



AC 



2a 2 
45m 4 



(E* - B 2 ) 2 + 7(E • B) 



(73) 



In analogy with what we did at finite temperature in the previous section, we now can 
obtain the correction to the Casimir energy from 



AE = — J d 3 x{ AC ] 



(74) 



in lowest order perturbation theory. 



Let us first assume that this energy correction can be written as an integral over a 
local energy density AS. Its magnitude will be given by AS = —{AC). Using Wick's 
theorem in the form ( E 2 E 2 ) = (E 2 }(E 2 ) + 2( Ef ) 2 5ij with the expectation values from 



(|6S|) and ( |70[) we then obtain 



AS 



2a l 



15m 4 



5(E' x y + 5(E 



2 \2 



{E 2 X ){E 2 



a 2 vr 4 



2 7 3 3 5m 4 L 8 



11 

225 



+ 9F 2 (6) 



(75) 



Because of the small, numerical coefficient in front and the condition L ^> 1/m this 
correction is seen to be negligible for all practical purposes. However, its magnitude gets 
larger as one approaches one of the plates where it is seen to diverge like AS(z — > 0) = 
— a 2 7r 4 /1620m 4 z 8 . The integral of this energy density over the space between the plates, 
will thus also diverge. Such a result is meaningless and a forces us to question the physical 
content of this position-dependent energy density. 



If we instead first do the volume integral in (|74j) and afterwards regularize the result, 
we get a finite result for the Casimir energy correction. It results from terms with the 
typical behavior of 



f L dz(Ef(x) 
Jo 



3/ dz(E 2 z ( X ))(E 2 z ^)[ 
/o 

2 




dz ^2 m 3 n 3 [l + cos (2nzm/L)][l + cos (2irzn/L)] 



m,n=l 
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The constant in the integrand simply gives C 2 (-3) = (1/120) 2 while the integrals over the 
cosines vanishes because of the boundary conditions. The last integral, from the product 
of the two cosines, is only non-zero when m = n. But then its value will be multiplied 
by C( — 6) = and is thus also zero. So it is only the position- independent parts of the 
field fluctuations that contribute to the total energy. The net result for the integrated 
correction to the vacuum energy is therefore given by just the first term in (|75[), i.e. 



lk 2 f 4 

M = -2^5W < 76 > 

Dividing by the plate separation L to get an energy density, we see that it can also be 
obtained from the last term in the free energy (|24] ) by the substitution T — ► 1 /2L just as 
the leading term of the Casimir energy could. It represents a tiny, additional contribution 
to the attractive force between the plates. 

The lack of commutativity between integration and regularization we have just ob- 
served in the calculation of the vacuum energy, is already present in the case of free fields 
as first pointed out by Deutsch and Candela sp^ . It can be seen by integrating for ex- 
ample the electric field fluctuations (|65| ) and fl69|) over the separation between the plates, 
which gives the electric half of the Casimir energy ( |62| ) after regularization. Only the 
position-independent portion of the fluctuation contribute to the result. On the other 
hand, integrating first the regularized field fluctuation in (|7l|) obviously gives a diver- 
gent result. The regularization dampens out the effect of the high energy modes of the 
field which prevents us from making any reliable predictions near the plates. As a conse- 
quence, one is not allowed to integrate a regularized result over the full volume between 



the plates. As shown by Scharnhorst 1 4C ] and Barton |4l[| the same vacuum fluctuations 
will also modify the speed of light in the space between the plates. But again one finds 
that the position-dependence does not contribute to the net effect. 

An unambiguous result for the energy density and pressure between the plates can 
be obtained from the energy-momentum tensors (^?]) and in the same way as was 
done in the previous section for the photon gas at finite temperature. In the expectation 
value of the free Maxwell part the position-dependent terms again cancel out and gives 
the traceless result 

(T$) = 7^(^-4^) (77) 

Here the four-vector = (0; 0, 0, 1) represents the normal n to the plates and we are 
using a metric rj^ u with positive signature. The energy density £ = ( Tqq ) is obviously 



consistent with the integrated result (32) while (T zz ) gives the Casimir force. Similarly, 
the contribution ( TjjjJ 1 ) from the interactions is fairly straightforward to calculate using 
the same field correlators as above. However, the evaluation of A( T^l ) in (^5|) is now much 
more difficult than at finite temperature. One must use the full field propagator discussed 
in Appendix C where the transverse and longitudinal components behave differently. But 
the main reason for the complications is the lack of momentum conservation in the z- 
direction which results in highly divergent, oscillating double sums. The evaluation of 
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these will require a separate investigation. At this stage we can only say that the energy- 
momentum tensor for the interacting system must satisfy the trace condition Tjf = —4AC 
which follows from (|37| ) when taken to Minkowski space. The expectation value of the 
right hand side is here the energy density AS in (|75|), 



4 Discussion and conclusion 



Effective field theory is a very powerful framework for calculating higher order quantum 
effects of interacting fields below characteristic energy scale. In our case of low-energy 
photons this is set by the electron mass. The Stefan-Boltzmann energy increases with 
temperature like T 4 and one would expect for dimensional reasons that the first radiative 
correction would go like T 6 /m 2 . We find that this term is in fact absent at low temper- 
atures since the Uehling interaction can be transformed away in the absence of matter. 
The first correction comes in at next order and is due to the Euler-Heisenberg effective 
interaction. It gives a contribution proportional to T 8 /m 4 . In the effective theory it re- 
sults essentially from a simple 1-loop calculation, while in full QED it would result as the 
low-energy limit of a 3-loop calculation. The absence of the T 6 term in the energy den- 
sity is also observed for the interacting pion gas calculated in the chiral limit of massless 
fields [42 1 . If there is any connection here, the reason is not clear. 



We have also calculated the finite-temperature polarization tensor of the interacting 
photons. In the effective field theory approach it is obtained from a 1-loop calculation 
while in full QED it would come from a more difficult 2-loop calculation. The resulting 
permittivity and susceptibility gives an effective light velocity which is less than one. At 
normal temperatures one can safely forget the small correction, but it will increase for 
extremely high temperatures which one has in the early universe. If it has any practical 
consequences in this context, does not seem to have been investigated. But it is clear that 
the propagation of light over cosmological distances is affected by the finite-temperature 
quantum fluctuations in the vacuum. 

Much of these finite temperature results can be applied directly to the Casimir effect 
induced by two parallel plates held at zero temperature. The absence of the T 6 term in 
the free energy implies that the first quantum correction to the Casimir energy density is 
caused by the Euler-Heisenberg interaction and varies like l/m 4 L 8 with the plate separa- 
tion. If the boundary conditions of the electromagnetic field on the plates could be taken 
to be periodic, one could in fact obtain the Casimir energy simply by letting T — ► 1/L in 
the finite-temperature results. The field fluctuations ( E 2 ) and ( B 2 ) would be constant 
between the plates and there would be no divergences as one gets near one of the plates. 
But the physical boundary conditions we use forces some of the components of the fields 
to be exactly zero at the plates. This classical constraint on the fields induces large quan- 
tum fluctuations which in fact diverge close to the boundary. The corresponding electric 
and magnetic energy densities diverge also and thus also their integrated contributions 
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to the total energy. For the free field, the divergences cancel so that one is left with a 
finite total Casimir energy. But for interacting photons, the energy density is in general 
divergent near the plates. In spite of this, we have obtained a finite Casimir energy for the 
interacting system by integrating the energy density first and then regularizing the result. 
We see that he processes of integration and regularization do not in general commute with 
each other. 

In order to analyze this apparently paradoxical situation a bit more carefully, let us 
look at the divergent sum 

oo 

S(9) = J2 sm 29n (78) 

n=l 

which enter the calculation of the basic correlator (|65|). The field fluctuations are given 
by its third derivative resulting in the function F{9) in (|66| ) when we use zeta-function 
regularization. When 9 = or 9 = n each term in the sum is zero while the regularized 
sum is infinite. So it cannot be used on the plates. A more physical understanding of 
the sum is obtained by evaluating it instead with an exponential cutoff as a regulator. 
Defining the convergent sum 

oo 

S( e ,9) = 2 e- en sm29n (79) 

n=l 

we will then have S(9) = lim e ^o S(e, 9). The cutoff should go to zero in the limit where the 
electron mass becomes very large since it is m which sets the lower scale for the frequencies 
which should be dampened out of the problem. Since the regulated sum is a geometric 
series, it is easily found to be 



e" e sin 29 
1 - 2e" £ cos 29 f < 

Expanding this result in powers of e one finds 

1 „ 1 cos6> 

- COt 9 ; 

2 8 sin 



S(e,9) = -cot 9- - e 2 + 0(e 4 ) (81) 



We see that as long as s'm9 > one can take the limit e — > and the first term will 
represent the fluctuation sum S(9) as we have done in the text. However, near the plates 
where sin9 — > the higher order terms in ( |8l| ) will become more and more important 
and the first term will no longer alone represent the fluctuations. The higher order terms 
with finite e must then be kept and represent in some way counterterms which must be 
introduced when one approaches one of the plates. 

In particular, this applies when we want to calculate the total energy between the 
plates. Then we are faced with the integral 




e~ e sm29 



2e" e cos 29 + e~ 2t 



J o 
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On the right-hand side we must take the derivative while e is kept finite. The integral 
gives then zero at the upper and lower limits. This is consistent with what we found in 
the text that the integral of the position-dependent energy density is in fact zero and not 
infinite. But it obtains only when we keep all the higher order terms in (|8l]) since we need 
the fluctuations on the boundary when we do the integral. 

In this way we get a finite result for the integrated Casimir energy also when photon 
interactions are included. All the position-dependence from the correlators integrate out 
to zero and we are left with a final expression which can be obtained simply by letting 
T — > 1 /2L in the corresponding result for the finite-temperature free energy of the photon 
gas. 

We want to thank E. Braaten, L.S. Brown and K. Scharnhorst for useful discussions, 
suggestions and clarifying comments. We thank NORDITA for their generous support 
and hospitality during the completion of this work. 
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Appendix A 

In order to illustrate how equations of motions can be used in effective theories, let us 
consider a scalar field <f> coupled to fermions if) and described by the effective Lagrangian 

C = ~<f>U<f> - gc^lf) + ^ ( D 0)2 

Since we want to simulate effective Maxwell theory in the text, the scalar field is taken to 
be massless. The last interaction term where a is some dimensionless coupling constant 

□ 2 

• 

Figure 8: Derivative correction to the boson propagator. 

and M is a heavy mass, is meant to be the corresponding Uehling interaction. It gives a 
correction to the scalar propagator shown in Fig.||. 

The Yukawa coupling g gives rise to a fermion-fermion interaction in lowest order g 2 
due to an exchange of a massless 0-particle. Because of the above Uehling term, there 
will be corrections to this Yukawa interaction. To lowest order in a it results from the 
Feynman diagram in Figj9| which gives the effective vertex 

which is a local four-fermion interaction. 
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Figure 9: Propagator correction to the four-fermion interaction. 



This result can also be obtained by using the equation of motion □(/> = g^t/j for the 
scalar field which results from the dimension D = 4 part of the Lagrangian. Putting now 
this back into £, we obtain 

.1 - a - 

£ = --(f>U(j) - gcjnpip + —g^Ucf) 

when we keep only interactions of dimensions D = 6 or less. In this way we have now 
generated a new derivative Yukawa coupling between the fermion and the scalar field. 
When combined with the original Yukawa coupling, we obviously get the same, effective 
four-fermion interaction as above now due to exchange of massless (f) particle with the 
standard 1/k 2 propagator. Instead of using the equation of motion, we can perform the 
field redefinition 

in the Lagrangian C. Keeping only terms to order 1/M 2 , we see that we then recover the 
modified Lagrangian £'. 

When there is no fermionic matter present, the Uehling interaction can thus be trans- 
formed away. Hence, it should also then give no contribution to the vacuum energy. This 
can also be seen by considering the lowest order Feynman diagram in Fig.[l(] resulting from 
the original Lagrangian C It gives a contribution 



M 2 J (2vr) 4 k 2 

when we use dimensional regularization. Had we instead used a momentum cutoff A, it 



Figure 10: Lowest correction to the vacuum energy. 

would have been oc A 6 and would be cancelled by an appropriate counterterm. Also at 
non-zero temperature this contribution to the free energy can be set equal to zero. 

If we denote the scalar propagator by D(x — y), the same diagram in coordinate space 
would have given 

A£ = ^yW nlD{X - y) = W 6 (0) 
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which obviously is the same as the above result. But the point now is that if the scalar 
field is confined by some bounding surfaces, its propagator will instead be changed into 
D(x,y). The Casimir energy density due to the the same diagram is then given by 



A£=lim-^a 2 x D(x,y) 



x-+y M 2 



But even the confined propagator satisfy the basic equation D x D(x, y) = S(x — y) so that 
the contribution to the vacuum energy density is the same as for an unconfined field and 
obtained above. The Uehling contribution to the Casimir energy is therefore zero. 



Appendix B 

We will here present the electromagnetic field correlators at finite temperature and de- 
rive some of the summation-integrals which appear in evaluating the relevant Feynman 
diagrams. 



From the basic, free field correlator ( |13D we can immediately write down the corre- 
sponding electric propagators, 

( Bk(x )E j {y) ) = ^ e^-y) (82) 

the magnetic propagators, 

( B^x) B 3 {y) ) = ^ z lHx - V) (83) 

and the mixed propagator 

<^(x)5 i (y)>=5fM^e*(»-») (84) 

Notice that this latter one is zero when x — > y so that the expectation value ( E • B ) = 
consistent with parity invariance. 

The electric field fluctuation ( E 2 ) now follows from the propagator fl82| ) by taking 
the coincidence limit x — > y. We then recover the result Jl?] ) in the main text. It can be 
evaluated from the basic integrals 



and 



r f.2 

hi = % i, , ,2 (86) 
% k 2 + ki 
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both of which can be obtained from the master integral (61) in the text. From rotational 
invariance we have that Lj = 5ijI kk /3 so we need only I kk . With dimensional regulariza- 
tion it follows that I44 + I kk = and one is left with only the integral -Z44. Then it follows 
that 

~ ^ {27TnT) 2 00^ (2 ^ T)rf r(1 _ d/2) 

44 ^7 (2^ k2 + (27mr)2 n ^ (4vr)^ T(l) 

There are no singular poles in the limit d — ► 3 where we find 

00 

I 44 = 2T 4 vr 3 / 2 r(-l/2) £ n 3 = -vr 2 T 4 /30 

n=l 

using zeta-function regularization for the divergent Matsubara sum with £(— 3) = 1/120. 
Together with the corresponding result for Jy we then have the basic integral (|i~6|). 



The full propagators enter the calculation of the extra term A(T44 ) in (|35|), Let us 
consider the very first term which can be expanded to give 



/ A(E 2 (x)E 4 (y)> = 4/ d 4 y(E i (x)E j (y))(E i (x)E j (y))(E k (x)E k (y)) 
Jo Jo 

+ sf d 4 y(E t (x)E J (y))(E t (x)E k (y))(E J (x)E k (y)} (87) 



Now the last factor in the first term is just ( E 2 ) which we already know. Integrating over 
y makes the momenta in the two first propagators in the same term equal and we have 

^d*y(E l (x)E 3 (y))(E l (x)E 3 (y))(E k (x)E k (y))=^ k ( ^±^ j (E 2 ) 

Here we need several new integrals which again can be calculated with the use of ((61 
The first one gives 

r kf (2vrnT) 3 r(l/2) 2 , . 



when done like I44 above. Now the integral 



j4 = % (P + k\f m 



follows from the identity J44 = I44 — J4 which gives J4 = — 7r 2 T 4 /20. Similarly, the integral 

can be obtained from J4 = I kk — J which gives J = 7r 2 T 4 /12. Using these results in fl8S|), 
we then have 

f3 d 4 y(E t (x)E J (y))(E t (x)E J (y))(Ek(x)Ek(y)) = ~ 2 ^\ 2 (91) 
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for the first term in (|87|). 

The second term can now also be obtained using the same integrals, 

r/3 



<?y ( Ei (x)E j (y) > ( E t (x)E k (y) ) ( E j (x)E k (y) > (92) 

% (W + k 2 ) 2 '% p2+p2 " 2-33-52 ^ 

Combining the two parts in (p7|), we find 



p 1 n7r 4 T 8 

A(E 2 (x)E 4 ( y )) = -i^ (94) 



in agreement with the first term of (36) in the main text. The other terms can then be 
derived along the same lines. 



Appendix C 

From the quantization of the photon field between two parallel plates in Section 3, it is 
straightforward to construct the full propagators in the Coulomb gauge. For the electric 
field we find for the x-component 

2 fd 2 k T fdko i (k k y ) 2 + (k x k z ) 2 



sin{k z z) sm(k z z'y kT < XT -^- iko ^ (95) 



where K 2 = k 2 — k 2 -, — k 2 with k z = nir/L. The expression for the y-components has the 
same form with k x and k y interchanged in the numerator while the propagator for the 
z- component is 

2 fd 2 k T fdk ~ ik? 



T 



L J (2vr) 2 J 2vr ^ K 2 + ie 
cos(k z z) cosik.z'y^^-*'^- 1 ^-^ (96) 

The magnetic propagators have the same general structure. There are also non-diagonal 
propagators between different components of the electric and magnetic fields. 

When the two points in the propagators approach each other we get the field expec- 
tation values. The integral over ko can then be done by contour integration. For instance, 
for the ^-component of the electric field we find from (|95[) 



<^ (X)) = z/tSI £ ^sm 2 (M (97) 

' n=l 
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where now uj = yk^ + fcf. This is in agreement with the corresponding expression (|6J 
in the main text. 



Bordag et al |15(] have succeeded in deriving a compact expression for all these field 
propagators based on quantization in the Lorentz gauge d^A^ = 0. With the plates normal 
to the z-axis at the coordinate values z = ao and z = a%, they found 

(M*)MV)) = iD^,v) = \ J 0f^e iPa{X - y)a e trlX3 -^ (> % <''''*' " 
with d 3 p = dpodp\dp2 and p a x a = PoXq —p\X\ — P2^2- Here P^ v is the projection operator 

9^- V -fr for /z, 1/ ± 3. 
for = 3 or v = 3. 



and r(p) = yPo — Pi — Pi- ^he position of the plates enter also through the matrix hij 
with the inverse 

/ e -«r|ao— oi| _j 



2 sin r | ao — a\ 



\ 



.\ e — ir|oo— ai| 



In the following we will take ao = and a\ = L as in the main text. 

The field expectation values can now be obtained by taking different derivatives of 
the propagator and then let the field points coincide. In order to regularize the ensuing 
divergences, we take the fields in the two points (x a , z) and (y a , z+e) where we let x a — > y a 
and leave the regulator e infinitesimally small, but finite. For instance, when we calculate 
the fluctuations in the x-component of the electric field we encounter the expectation value 

J_ = (doA^Ax)^ Ymx a ^—^-{A l {x)A 1 {y)) 



2 J (2rfr r )i 3 e 



2 J (2tt) 3 r 

By the Euclidean transformation p$ — » ip^ under which T — » i'y this becomes 
If dp i dp 1 dp 2 pl I pi 



1-a e -Tk-«d ^-1^ e -7k+6-a 3 -| 



2 7 (2?r) 3 7 ^ 7 : 
where now < z < L. Introducing spherical integration variables we then obtain 

1 roc ,-,,3 

15^ 2 i 2sinh( 7 L) L 
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after integrating over angles. While the function J_ is needed for evaluation of (E%), (Ey) 
and (Bg), we have here also given the related function I + which will be used to calculate 
(Ef), (B%) and (By). The remaining integrals are given in ||45|1 , 



dx 



2 1 -»r(n)( H (v,(p + i)/2) 



sinh(x) 

where Ch(s,x) is the Hurwitz zeta-function 



( H (s,x) = J2 7 — : — v 



Thus, we have the results 
1 



I± 



[Or (4, z/L) ± Crr(4, e/2L) ± Ch(4, -e/2L) + Ch (4, 1 - z/L)] 



40vr 2 L 4 

In the limit e — ► 0, it follows from the above definition that 



( H ( s ,e) = - + ((s) + 0(e) 



when s > 1. We then have 



I± 



1 



QO/~ 4 ~-4 

( H (4,z/L) + C H (4,l- Z /L)±-^±- 



For integer s = k the remaining zeta-functions can now be expressed in terms of the 
digamma function by the relation 



( H (k,x) 



-l) k d k ~ l 



(k- l)\dx k ~ 

In addition, using the identity ip(x) — ip(l — x) = — 7rcot(7ra;), it follows that 

, , / 1 -i \ Id 3 , , vr 4 / 3 2 

Cff(4, *) + 0(4, 1 - x) = kcotM] = T I - 



With 9 = irz/L as in the main text we can then write the result as 



4 7T 2 
i± = ±— — + 



"5vr 2 e 4 120L 4 



(98) 



where the function F(9) is defined as in fl67|). This gives (OqAiOqAi) = I_. Similarly we 
obtain (c^Ag^i^o) = ^OqAiOiAq) = ^— • Combining such terms, we then have 



(E 2 X ) = (OqA^A, + 28^3^ + c\ A)> = ^- 



7T 



48L 4 



1 

15 



+ F(9) 



28 



when we discard the divergent, distance-independent term. Both (Ey) and (B^) can be 
evaluated in terms of For example, in the calculation of (B%) we need the expectation 
value 

(d 1 A 2 d 1 A 2 ) = lim (A 2 (x)A 2 (y)) 
i"^" dxi dm 



2 J (2vr) 3 7 ^ 7 2 

In this way we find that (E^.) = (Ey) = — (B%) which agrees with what is obtained in the 
Coulomb gauge. Exactly along the same lines we also derive (El) = — (B%) = —(By) in 
terms of I+, i.e. 



(El) = lu ?2 



2 r 48L 4 



when we again drop the divergent piece which represents the field fluctuation in empty 
space without any disturbing plates. 
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